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Abstract 

We study the dynamics of a finite chain of diffusively coupled Lorenz oscillators 
with periodic boundary conditions. Such rings possess infinitely many fixed states, 
some of which are observed to be stable. It is shown that there exists a stable fixed 
state in arbitrarily large rings for a fixed coupling strength. This suggests that coherent 
behavior in networks of diffusively coupled systems may appear at a coupling strength 
that is independent of the size of the network. 

1 Introduction 

Lattices of coupled dynamical systems have been studied in many different contexts: as dis- 
crete versions of partial differential equations of evolution type |l(| , 1 22 , models of neuronal 



networks g , |18 and phase lock loops |4j , and in statistical mechanics |LJ] . The particular 
problem of synchronization and emergence of coherent behavior in lattices of diffusively 
coupled, chaotic, continuous time dynamical systems has received much attention recently 
due to its applications in neuroscience , Q , and chaotic synchronization [jl9| , Q . 

Analytical results in the literature suggest that the coupling strength necessary for the 
appearance of coherent behavior in networks of chaotic systems grow with the size of the 
network || , Q , jisj . For instance, estimates of Afraimovich and Lin j^] suggest that the 
coupling strength necessary to synchronize a lattice of generalized forced Duffing systems 
grows as a high power of the size of the lattice. Since stable, coherent states in lattices 
may frequently be viewed as instances of generalized synchronization Q, one might guess 
that the appearance of dynamically stable states of this type also requires that the coupling 
grows with the number of oscillators in the system. This would make it unlikely that such 
networks are of importance in nature, since the coupling strength in realistic systems cannot 
be varied by orders of magnitude. 

In contrast to these theoretical results, the numerical experiments of Huerta, et. al. [jl8| 
indicate that partial synchronization may appear in certain networks of chaotic systems at 
a nearly constant coupling strength, regardless of the size of the network. In particular, 
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Jig ] considered a system with about 10 coupled neurons. Each neuron was modeled by a 
system of three nonlinear, ordinary differential equations which, when uncoupled, underwent 
chaotic motion. They observed a variety of different stable, coherent motions including a 
number of states that exhibited roughly periodic patterns. It is the goal of this paper to 
show rigorously that such such stable, coherent states can occur in arbitrarily large systems 
of coupled chaotic oscillators at a coupling strength which is independent of the size of the 
system. 

Since the Lorenz system is one of the paradigms of chaotic flows, we have chosen a ring 
of diffusively coupled Lorenz equations for our investigations. The systems in the ring evolve 
according to the equations 

x\ = a(y l - Xi) + d x Ax t 
y'i = rx t -yi- XiZi + dyAyi (1) 
z'i = -@Zi + Xiin + d z Azi 

where Axi — Xi—i — 2xi + Xi+i with the index i taken modulo n is the discretized Laplacian 
with periodic boundary conditions. The constants a = 10, r — 27, and (3 — 8/3 are chosen 
so that each uncoupled system would be in the chaotic regime. Following the work of ]l8| , 
we will focus on the case d x ^ 0, d y = d z = 0. As shown in Q, such a ring is synchronized 
identically for sufficiently large values of d x . We will be concerned with the case when d x 
is sufficiently large for coherent spatial structures to emerge in the chain, but insufficient to 
synchronize the ring. 

The paper is organized as follows: Typical spatial structures observed in numerical 
experiments are described in section ^. In particular, we give examples of stable stationary 
states, traveling waves, and breathers. In section || a general framework for the study of fixed 
states in chains and rings of dynamical systems is discussed. We show that the existence 
of such states can be reduced to the study of periodic solutions of reversible dynamical 
systems on 1™. We then apply these results to two different cases, including the chains of 
coupled Lorenz oscillators (|l|) , to show that a large number of fixed states can be expected 
in such rings. The conditions under which these states are stable are investigated in section 
H We then show how Floquet theory can be used to reduce the stability question in lattices 
of arbitrarily many oscillators to the study of the spectrum of a matrix of fixed size (but 
depending on a parameter - the "quasi- momentum" .) Finally, in section ^| we implement 
this program to rigorously construct an example of a spatially periodic stable fixed state 
which remains stable for rings of arbitrary size. 

2 Numerical Experiments 

Numerical experiments were performed on a Sun SuperSparc using the numerical integration 
program XPP. Several different integration methods were used to check the validity of the 
results. These included fixed and variable step Runge-Kutta methods of varying order and 
the Gear method. 

At very low values of the coupling constant d x , the individual systems in the ring behaved 
as if they were uncoupled. As the coupling strength was increased the behavior of the 
ring became more coherent. Since the discrete Laplacian acts to synchronize neighboring 
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systems in the ring, these results agree with our expectations. As was emphasized in the 
introduction, the onset of the coherent behavior occurred at coupling strengths which seemed 
to be independent of the number of oscillators in the ring. Some of the typical structures 
observed in these numerical experiments are discussed below. 

2.1 Breathers 

If some of the systems in the ring are nearly stationary, while others are undergoing oscil- 
lations we say that the oscillating systems form a breather. Since we are considering rings 
of finite size, it is not expected that any single system in the ring is stationary, unless the 
entire ring is stationary. 

An example of a system with two breathers is shown in figure It is interesting to note 
that in these states each of the individual systems undergoes oscillations that are nearly 
periodic, and that there is a strong spatial dependence on the motion. Loosely speaking 
the motion of each individual system is trapped in a region close to one of the 'lobes" 
of the Lorenz attractor. No stable states in which all systems in the ring are trapped in 
the vicinity of only one "lobe" has been observed. Similar breathers have been studied in 
coupled map lattices in Q , while occurrences of nearly periodic behavior with pronounced 
spatial patterns in lattices of chaotically bursting neurons have been observed numerically 
by M.I. Rabinovich, et. al. 

Another interesting feature of breathers, which is also shared with other types of solutions 
of this network, is that frequently nonadjacent systems synchronize without synchronizing 
with the oscillators that lie between them. The synchrony is usually not exact, however 
this may be due to numerical errors in the calculations. Due to the many symmetries 
of the system this sometimes seems to be due to the stability of an invariant manifold 
corresponding to a partially synchronized state in the chain as discussed fl9[| . In other cases 
the state does not appear to be symmetric, and it is unclear what the mechanism behind 
this synchronization is. 

It is important to note that some of the qualitative features of the dynamics remain 
unaffected by the size, or the coupling strength in the ring. Figure ^ shows typical timeseries 
of the x variable at four different coupling values in a ring of 32 oscillators. Note the smooth 
envelope and regularity of the oscillations, as well as the agreement of the timcscales among 
the different timeseries. These general features remain unchanged as the number of systems 
in the ring is varied, as long as the coupling is below the synchronization threshold and 
above the value necessary for the appearance of coherent behavior. 

2.2 Stable Stationary States 

In numerical experiments with rings of 8 to 52 systems stable steady states were observed 
for a variety of initial conditions for d x > 10. Typical stable fixed states are shown in figure 

I 

For certain values of the coupling it appears that all initial conditions lead to one of the 
stationary states of the system. The basins of attraction of the stationary states seem to 
be intertwined in a complicated way in this shown in figure |j. The same is true for 

other stable states of the system. Small changes in the initial values of the parameters, and 
changes in the integrating methods can result in very different asymptotic behavior of the 
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system. This situation is similar to that of attractors with riddled basins of attraction , 
and was also observed in numerical simulations of networks of chaotically bursting neurons 

2.3 Traveling Pulses 

Traveling pulses were observed in the ring only when d x ^ and d y ^ 0. A typical traveling 
pulse is shown in figure B. The pulse oscillates as it propagates along the chain, and is thus 
a periodic, rather than a fixed state in a moving coordinate frame. This can be seen in the 
timeseries in figure 

3 A General Framework for the Study of Steady States 

System ([j]) is a special case of a lattice discrete or continuous time dynamical system 

K)' = f(W s ) (2) 

where j £ 1 D , Uj £ M. p and {uj} s = {ui : \i — j\ < s}. In the following discussion we will 
use the convention that in u\ the subscript i £ TP denotes the position in the lattice, while 
the superscript 1 < k < p denotes the component of the vector m. Such systems have been 
studied by many authors M , , Q as discrete versions of partial differential equations 
of evolution type. 

We will mainly be concerned with systems continuous in time with D = 1 and s = 1. 
These are simply chains of systems coupled to their nearest neighbors. The special case of 
rings, that is chains of finite size with periodic boundary condition will be the focus of our 
attention. A state in a ring of n systems corresponds to a state of period n in the spatial 
variable in an infinite chain. 

A steady state of a chain given in (|^) is determined by 

F({u j } s )=F(u j . 1 ,u j ,u j+1 ) = (3) 
F ({ujY) = F{u ] - ll u ]1 u ]+1 ) = Uj (4) 

in the continuous, respectively discrete case. We will study systems that satisfy the following 
condition: 

Condition 1 If we write F : R 3p — > W as F(x,r),Q = (F 1 (x,V,0, F 2 (x,r),C), 

- s J 1,3 

By the Implicit Function Theorem in the case of continuous time if det §tj-(X; V> C) 7^ 

L ^ J i.j 

for all values of x, Vi C then either there exists a point (a, b, c) £ R 3p such that F(a, b, c) = 
and therefore a function G(r],x) such that F(x, r), G(r],x)) — 0, or F(x, f], C) = has no 
solutions and there are no fixed points of (|^). If such a function G exists and we define 
Uj + i = G(uj, Uj-i), Xj = Uj and yj = Uj-\ this leads to the following dynamical system on 

x j+ i=G(x j ,y j ) (5) 
Vj+i = x 3 
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The steady states of (||) are given by the x-coordinates of the orbits of (ph. In general, the 
function G may not be unique, in which case more than one system of the form (^|) is needed 
to determine all the fixed points of the chain. An equivalent argument holds in the case of 
discrete time. 

Since the function G can assume any form, not much can be said about such systems in 
general. The subclass of chains of systems with the following type of coupling is easier to 
analyze. 

Definition 3.1 A nearest neighbor coupling of a chain of systems is said to be symmetric 
if F(uj-i,Uj, Uj+i) = F(uj+i, Uj, Uj-i). 

This simply means that a system in the chain is coupled to its left and right neighbor in the 
same way. For instance the couplings Auj = Uj-% — 2uj + Uj+i and i&Uj = Uj-%Uj+x + Uj 
are of this type. The stationary states of symmetrically coupled chains are related to the 
following class of dynamical systems: 

Definition 3.2 Given a diffeomorphism <£> : R 2p — ► M 2f> and an involution R : R 2p — > R 2p 
such that the dimension of the fixed point set of R is p we say that $ is R-reversible if 

i?o$ = $- 1 oi? (6) 

The dynamical system defined by Xi + \ — $(2^) is also said to be R-reversible. 

The relation between symmetrically coupled chains and reversible systems is given in the 
following 

Theorem 3.3 The fixed states of a symmetrically coupled chain of systems satisfying con- 
dition 1 correspond to the orbits of a dynamical system of the form (^) which is R-reversible 
and volume preserving. For such a system R(x, y) = [y, x), with x, y £ W. 

Proof: The arguments for the continuous and discrete case are virtually identical, so only 
the first will be considered. Equation (^|) together with the assumption that the coupling is 
symmetric implies that 

F(uj- 1 ,v,j,iij +1 ) = F(u j+li Uj,Uj-x) = (7) 
Therefore this leads to the recursion equations 

Uj+i = G(uj,Uj-i) Uj-i = G(uj,Uj+x) (8) 
and the two dynamical systems 

x J+1 = G(x v y 3 ) Zi _ x = G(z v w 3 ) 

Vj+i = 1 = z i 

where and Wj = Uj+i- The action of these systems is shown 

schematically in figure H. 

Let R(x,y) = {y,x) and $(x,y) = (G(x,y),x) so that the dynamical systems in (S) is 
generated by <f>. By the definition of this diffeomorphism 

Ro $ o Ro $>(x h yi) = Ro $(y i+l ,x i+l ) = Ro<^(z l ,w l ) 

= R(zi-i,Wi-i) = R(yi,Xi) = (xi,yi) 
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Since R is an involution with Fix(i?) = {(x,y) : x = y} the diffeomorphism $ and the 
dynamical system it induces on the plane are i?-reversible. 
To prove that the map <J> is volume preserving notice that 



D<5> 



D\G D 2 G 
I 



(11) 



By the definition of the function G we know that F(x, Tj, G(rj, x)) — so that differentiating 
with respect to x leads to 



DxF + D 3 FD 2 G = 



(12) 



D3F is an invertible matrix by Condition 1, so that D 2 G = —DiF(Dj,F)^ 1 . By assumption 
the system is symmetric so that D\F = D^F at all points in M. 3p and hence D 2 G = —I. 
Since 



" I 




£»iG -I 




/ 0" 


-I 




I 




-L>iG / 



and 



det 



/ 
-I 



= det 



/ 

-D X G I 



= 1 



(13) 



(14) 



it follows that det — 1 and the diffeomorphism $ is volume preserving. <0> 

i?,-reversible systems have many properties that facilitate their study. We will make use 
of several of these summarized in the following proposition adopted from [fl2| 



Proposition 3.4 Let R and $ be as in Definition (3.i). 

a. If p E Fix(R) and $ fc (p) £ Fix(R), then Q 2k (p) = p. Such periodic points will be 
referred to as symmetric periodic points. 

b. LetpE Fix(R) be a fixed point 0/$, then R(W u {p)) = W s (p) andR(W s (p)) = W u (p) 
so that if q £ W w (p)n Fix(R) then q is a homoclinic point. 

c. If p £ Fix(R) is a fixed point 0/$ such that W u (p) intersects Fix(R) transversally at a 
point q, then there exist infinitely many symmetric periodic points in any neighborhood 
of p. 



4 Steady States In the Case of Discrete Laplacian Cou- 
pling 

For a chain of systems of the form u' = f(u) coupled through a discrete Laplacian Auj = 
(uj-i — 2uj + Uj+i) the equation (||) takes the form 

( Uj y = f(uj) + diuj-t - 2uj + u j+ i) (15) 



G 



and the dynamical system (||) takes the form 



1 c 

x l+1 = --f(xi) - y, t + {2 + -)xi (16) 

Vi+i = Xi 

where c = in the case of continuous time and c — 1 in the case of discrete time. The 
following proposition follows immediately from these definitions 



Proposition 4.1 The fixed points of ( jiq j are of the form xi = yt = uo where uq is any 
solution of the equation f(u) — in the continuous time and f(u) = u in the discrete 
time case. Thus all fixed points of ftlb] ) lie on Fix(R) — {(x,y) : x — y} and are in 1-1 
correspondence with the steady states of an uncoupled system in the chain. 

The fixed points ofjq) correspond to states of the chain that are constant in space and 
time and proposition |4.l| shows that the only such states {uj}j e % in the case of a discrete 
Laplacian coupling are given by Uj = uq for all j 6 Z where uq is a fixed point of an 
uncoupled system of the chain v! = /(it). Since the fixed points of Jl6| ) are [uq,uq) € 
Fix(i?) this allows us to make use of proposition 3.4.c[] 



Moreover we can apply Theorem 3^2 directly in this case to conclude that the map 
defining the fixed points is irreversible and volume preserving. 

Remark: It is interesting to note that if the fixed states of a chain of one dimensional 
systems are determined by a dynamical system of the form ( |l6| ) on M 2 and (xq, yo) is a fixed 
point of ( ft6| ) then the discriminant of the characteristic polynomial of DQ(xq, yo) is 



{-- d f'{x a ) + 2 + c/dy -4. (17) 

Since det D&(xq, yo) = 1 we can conclude the following: 

a. In the case of discrete time the eigenvalues of D&(xq, yo) are on the unit circle when 
f'{xo) > 1 and d > (f'(x Q ) - l)/4) or f(x ) < 1 and d < (f'(x ) - l)/4). 

b. In the case of continuous time, the eigenvalues of _D$(xo, yo) are on the unit circle for 
sufficiently large positive d when f'(xo) > and for sufficiently large negative d when 
f'(xo) < 0, 

By a theorem of Birkhoff (see for example p0|) since $ is area preserving, the dynamical 



system (16) will generically have infinitely many periodic orbits in any neighborhood of 
(xo,ya). Any such periodic orbit will correspond to a fixed state of ( [l5| ) which is periodic 
in space, i.e. in the variable j. 

The condition det D§(xo, yo) = 1 is not sufficient to guarantee that D&(xo, yo) will have 
eigenvalues on the unit circle if ( |l6| ) is not a dynamical system on the plane. Additional 

1 The chains under consideration are a special case of chains of the form 

= f( u j) + g( u j-i> u j, u j+i) 

where g(uj-i , u 3 ;, Uj+i) is a symmetric coupling vanishing on the linear submanifold of R 3p defined by 
Uj — i = uj = Uj+i = 0. Such systems can be analyzed using the approach of this section 
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information needs to be considered to conclude the existence of such periodic orbits in this 
case. 



Example 1: Consider a chain of tent maps coupled through a discrete Laplacian de- 
scribed by the equation u'j = f(uj) + Auj with 

I — 2uj + 2 11 Uj > 1/2. 

This is a discrete dynamical system with fixed states determined by a homeomorphism of 
the plane given in ([l6]) with c = 1 . 



The only fixed points of (|16[ ) are (0,0) and (2/3, 2/3) by proposition 4J. Since in this 
case the dynamical system (|16|) is linear in each of the two halves of the plane x < 1/2 and 
x > 1/2 the analysis of the dynamics around the fixed points is straightforward. A direct 
calculation shows that in the region x < 1/2 there is a family of invariant ellipses whose 
major axis lies on the diagonal D = {(x,y) G M. 2 : x — y} and on which the action of (16) 
is a rotation. As d is increased, these ellipses become more eccentric. 

The fixed states of the chain corresponding to these orbits are not dynamically stable. 
If u = {iij}j£z is a fixed state in the chain such that uj < 1/2 for all j then around this 
state the dynamics of the chain are described by: 

= 2uj + dAuj (19) 

The spectrum of the operator 2 + dA is the interval [2 — 4d, 2] and so u cannot be stable, 
and would not be observed in numerical experiments. 

Besides these invariant ellipses, this system will also typically have infinitely many pe- 
riodic orbits, and will exhibit complicated behavior, as the following argument shows. A 
simple calculation shows that the point p = (2/3,2/3) is hyperbolic for all values d > 0. 
Since the dynamics is piecewise linear the stable and unstable manifold can be calculated 
explicitly. W u (p) will consist of a ray contained in the half plane x > 1/2 and a second more 
complicated part which is constructed as follows. The first section of W u (p) is a line A, as 
shown in figure ^. The part of A contained in the half plane x < 1/2 is rotated along the 
invariant ellipses around the origin to the line A' under forward iteration. The subsequent 
images of A' are lines that are rotated further. Therefore f n °(A') intersects the diagonal 
D for some no. Whenever the angle of this intersection is not a right angle W u (p) will 
intersect W s (p) transversely since by proposition [TJ.b the stable manifold is the reflection 
of the unstable manifold through the diagonal W u (p). This shows that complicated dynam- 
ics of ( |l(j| ) can be expected and implies the existence of infinitely many periodic as well as 
spatially chaotic states in the chain of tent maps. 

Example 2: Next we consider the chain of diffusively coupled Lorenz system (P with 
d y = d z — 0. This system does not satisfy condition 1, however after setting the right -hand 
side of equations ([!]) to and using the second and third equation in (|l]) to eliminate the 
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variables y.; and Zi the following equation for a fixed state is obtained 



a(xi 



f3rxi 



P+(xi) 2 



y> 



f3rxi 



2,Xi 



= 



P+(x<) 2 
~/3+(^) 2 ' 



(20) 
(21) 
(22) 



The function on the left hand side of equation ( |20| ) satisfies condition 1 , and we can proceed 
as in the previous example. 

Equation (J20|) defines the following dynamical system on the plane 



Xi+\ 






)" 




Vi+i_ 




Vi. 







(23) 



Orbits of period n of this system are in 1-1 correspondence with the steady states in the ring 
of n Lorenz systems. The two fixed points of ( p3| ) are (0,0) and (±-^//3(r — 1), ±\J P{r — 1)) = 
(±8.326, ±8.326) which both lie on lie on Fix(i?). 

The following definition from |l2| will be used in the remainder of the argument. 

Definition 4.2 A compact region H in the plane is called overflowing in the y- direction for 
a function F if the image of any point (xo, yo) € int(H) lies strictly above the line y = yo. 

Notice that any point in the interior of H either leaves H or is asymptotic to a periodic 
point on the boundary of H. 



Theorem 4.3 The map (\23() has a homoclinic point and infinitely many periodic points for 
all0<d< 20. 



Proof: Let p denote the fixed point (— yj [3{r — 1), —y//3(r — 1)) of (|23]), G the ray 
= —y/f3(r — 1) with x > —yjHr — 1), and D = {(x,y) : x = y} the diagonal in R 2 . The 



y 

image of the triangular region W bounded by D and G lies between the two cubics 

Pry 



F(D) = {(x,y):x=-{y 



P + y 



:) + y} 



F(G)={(x,y):x = -(y 



Pry 
P + y' 



) + 2y+ y/P(r-l)} 



(24) 
(25) 



Consider the region F(W) n W = H depicted in figure |?]. Let (xq, yo) G mt(H) and consider 
the vertical line segment I passing through (xc>>J/o) and connecting G and D. The image of 
I is a horizontal line segment connecting F(G) and F(D) contained in the line y = xq. Any 
point of this line segment lies above the line y — yo. Therefore H is overflowing in the y 
direction. 

Since mt(F(H)— int(iJ)) and 'mt(H) lie on opposite sides of D, any point in H must 
have an iterate which either lies on D or crosses D. A direct calculation shows that one 
branch of W M (p) enters H. Since H is overflowing in the y direction, W u (p) must either 
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cross D or be asymptotic to the fixed point (0, 0) creating a saddle-saddle connection. The 
second possibility can be excluded as follows: 

The eigenvector of the linearization of F around (0,0) corresponding to the stable direc- 
tion is Vl = (1, - d _ 130+ 2V65"-rf )- The vector V2 = C 1 ' is tan S ent t0 F ( D ) at (0,0). 
A direct computation shows that for < d < 20 the vector vi points to the left of V2 
and so the tangent to W s (0,0) at (0,0) does not point into H. This situation is depicted 
schematically in figure || and excludes the possibility that W u (p) is asymptotic to (0,0). 

By proposition 3.4.b the fact that W u (p) meets D implies the existence of a homoclinic 
point. 

By proposition 3.4. c in the case of irreversible systems an infinite number of periodic 
points exists whenever W u (p) crosses D transversely at a homoclinic point. If this point of 
intersection is not transverse then a transverse intersection of W u (p) and D can be produced 
nearby by an argument given in [[[21 p. 261]. 

Numerical investigations suggest that the theorem remains true for arbitrarily large 
values of d. 

This theorem shows that rings of Lorenz systems coupled through the first variable can 
be expected to have many periodic stationary states. Numerical computations of W u (p) 
suggest that even more complicated stationary states can be expected (see figure ^J) . However 
it is unclear whether any of them are stable. Numerical investigations show that not only 
are they stable, but their basins of attraction occupy a large portion of phase space. We 
next address the problem of stability of these fixed states. 



5 Stability of Steady States in a Ring 

The study of stable fixed states in a chain of oscillators is not new. They are described in 
fill , H as instances of oscillator death in a chain of neural oscillators, in Q] as synchronous 
states of phase lock loops, and conditions for the stability of complex stationary states 
in general lattices are given in ||. The situation presented here is different in that an 
explicit periodic stationary state is analyzed in a case where the methods from the theory 
of parabolic partial differential equations used in Jl3| are not applicable, and the conditions 
proposed in || cannot be verified. In addition, since the coupling strength is neither very 
large nor very small, there is no obvious perturbative approach to the problem. 

Even if some of the fixed states in a ring of n systems are stable it can not be concluded 
that stable equilibria exist in arbitrarily long rings of such systems. If a stable fixed state 
is viewed as a special case of synchronization, one might guess from results in Afraimovich 
and Lin || that the coupling necessary to have stable steady states in longer rings increases 
as a power of the size of the ring. It has already been argued that this is not realistic from 
a physical viewpoint, and we will show below that it is not the case. 

In the following sections we will concentrate on the specific example of a ring of Lorenz 
systems coupled in the x variable discussed in example 2 of the last section. However the 
ideas presented can be applied to any chain of symmetrically coupled systems for which 
condition 1 holds. 

The linearization of each of the Lorenz equations around a point Vj = (xi,yi,Zi) when 
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(27) 



where the matrix T is defined as: 



r = 



d 





(28) 



It is in general not possible to compute the eigenvalues of this matrix analytically. As 
mentioned above, we are specifically interested in coherent states in networks with arbitrarily 
many oscillators. Given a periodic stationary state in a chain of n oscillators, an easy way 
to obtain a periodic stationary state in arbitrarily long chains is to repeat this state. In 
particular, if v is a fixed state in a ring of n oscillators then repeating these values K times 
to obtain v^ K > = {v, v, . . . , v} produces a steady state in a chain of Kn oscillators. We refer 
to such a state as the K -multiple of the steady state v. An example of a stationary state and 
its 4-multiple is given in figure [l^. We next derive conditions under which all if-multiples 
of a stable steady state are themselves stable and demonstrate that these conditions hold 
in a particular case. 

The linearization around the steady state leads to the following stability matrix: 



where 
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and E is an 3n x 3n matrix of the form 
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In general it is impossible to compute the spectrum of this 3Kn x 3Kn matrix. However 
in the present instance Bloch wave theory, provides a way of expressing the eigenfunctions 
of a larger system with a periodic dependency as the eigenfunctions of a smaller system. As 
usual, this reduction is accompanied by the introduction of an additional parameter into the 
equations. Here the Bloch wave approach will be used to reduce the problem of computing 
the eigenvalues of the 3Kn x 3Kn matrix D^k) F for arbitrary K to the computation of 
the eigenvalues of an 3n x 3n matrix dependent on a parameter t. 

Let e = V(j))f^i = (*(1), ri(l), *(Kn), $(Kn), n(Kn)) where 



= ex P 



Hj) =exp (^-j ) §,C?) 



T](j) = cxp 



V nK 
f 2iriq 
{ nK 



(29) 



3 Vq(j) 



for 1 < q < Kn and ?7 g are assumed to be n-periodic. We will show that the vector 

e is an eigenvector of D^k) whenever the vector 



i = (#(i),$(i),r?(i)) 



(30) 



is an eigenvector of a particular 3n x 3n matrix derived from D$( K )F. A direct calculation 
leads to: 

2iriq ■ 2iriq ~ 1-niq ~ 

[D^ K )F] e (3j) =e~^ J (de~^ ^ q (j + 1 mod n) + de~~^^ q (j — 1 mod n) — 
2*,(j) + «i + + M3)f}gU)) 

[D^(K)F] e (3j + 1) =e^'(u 2 (j)* g (i) + u*0')*«0') + u 2 {j)f) q [j)) 
[D^ K) F]e(3j + 2) =e^ j (v 3 (j)^ g (j)+w 3 {j)^ q (j)+u 3 (j)fj q {j)) 
for 1 < j < ifn. Here Vi(j),Wi(j),Ui(j) are given by equation (Eq) as 



(31) 



Vl(j) Wl(j) 

W2(j) W 2 (j) 142 (j) 

W3(j) w 3(j) U3(j) 



— CT <7 
^ mod n 1 

i/j mod n mod n 





mod n 

-b 



and are thus n-periodic. 

We want to check when e is an eigenvector of D^ K )F and so we set D^( K )Fe. = Ae. 
Using the expressions obtained in (|3l| ) we get 



exp 



/ 2itiq 



nK 



2-Kiq ~ 2-rviq ~ 

de~^^ ? (j + lmod n) + de ~^ q (j — 1 mod n) 



2*,(j) + wi(j)* 9 C?') + «»i0')*g0') + u i{j)n q ij), V2(J)%U) + Mj)^q(j) + M.i)v q {j), 



«3(j)*«(j) + W 3 (j)$q{j) + U 3 (j)fj q (j) 



Kr. 



J = l 



A exp 



(32) 
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Since all the functions on the right hand side of ( |32| ) are n-periodic we can express 
equation (p2n as T>{g)\ = Ai where i is defined in equation (pfl) and 



V(q) = 



D*J - 2T E(q)+ 
E(q)- D*J-2T E(q)- 







E(q)+ 
with the matrices E(q) ± defined as: 







de ± -^ L 



(33) 








(34) 



This implies that the eigenvalues of D^(k)F are the same as the eigenvalues of the 3n x 3n 
matrices T>(q) for 1 < q < Kn. The 3n x 3n matrix T>(q) is easier to handle than the 
3Kn x 3Kn matrix D^(k)F. 

To conclude that any if -multiple of a stable fixed state vo of a ring of n oscillator is 
stable it is sufficient to show that the eigenvalues of the matrix V(q) have negative real part 
for all values of q and K. To simplify the calculations we can replace the argument q G Z 
by a continuous parameter t G M by replacing the matrices E(q) :iz with the matrices 





~de ±lt 








Ef = 























in the expression for T>. If we can show that T>(t) has eigenvalues with negative real part 
for all £ G K then the same is true for T>(q) for any values of K and q in Z. 

The following lemma about the class of matrices of type (|33| ) will be used in the next 
section. 

Lemma 5.1 Let £ , (£) ± be as in the previous lemma. The characteristic polynomial of a 
Km x Km matrix 



D 



Mi E(t)+ 
E(t)- M 2 E(t)+ 







Eity 







E(t)+ 



E{t)~ Mk-i E{t)+ 
E(t)~ M K 



is of the form 



Kh 



(K-l)m 

p M (t,X)= +fc cos Kt)\ l + ^ 

i=(K-l)m+l 



(35) 



i=0 



A proof of this lemma is given in the appendix. 



13 



6 An Example 



This section gives a computer assisted proof that any if-multiple of a particular stable 
fixed state in a ring of four Lorenz system is stable. The proof involves the use of inter- 
val arithmetic and the following theorems about the convergence of the Newton-Raphson- 
Kantorovich method. The computations were performed using Mathematica's implementa- 
tion of interval arithmetic. 

Theorem 6.1 Let f{z) be a complex analytic function and assume f(zo)f'(zo) ^ for some 
z . Define h a = —f(zo)/f'(zo), the disc K = {z : \z — z \ < \hg\} and M = max^ \f"{z)\. 
If 2\1iq\M < \ f'(zo)\ then there is exactly one root of f in the closed disc Kq. 



Theorem 6.2 (Kantorovich's Convergence Theorem) Given a twice differentiate func- 
tion f : M" — * K" which is nonsingular at a point xo — (x\, x%, ■ ■ ■ , x„) let \Tik] — 
[jfi-fxo)]^ 1 ■ Let A,B,C be positive real numbers such that 



n n 

:^|r ifc |<A maxJ2\ T Mx )\<B C< — 



k=l 



k=l 



2AB 



(36) 



Define the region R = {xE M™ : max ! \x, - x° t \ < (AC)- X (1 - Vl - 2 ABC)}. If 

d 2 h 



max 



EE 

j=i k=i 



dxjdxi 



< C 



i = 1, . 



(37) 



then the equation f{x) — has a solution in R. 



The argument will proceed as follows: First a steady state of the ring is found using 
interval arithmetic. Next it is shown that the characteristic polynomial px>(\,t) of the 
matrix T>(q) in ([}3|) corresponding to this steady state has roots with negative real part for 
t = j. Finally it is proved that pp(A,£) does not have roots on the imaginary axis for any 
t £ R. Since the roots of a polynomial depend continuously on its coefficients this implies 
that the roots of p-p(A, t) cannot cross into right half of C. By the results of section || this 
implies that all iiT-multiplcs of the fixed state under consideration must be stable. 

In the following intervals of numbers will be denoted with overbars to distinguish them 
from real numbers, and make the notation less cumbersome. For instance A denotes an 



interval and the interval [3.16524,3.16533] is denoted 3.165. 

Numerical investigations with the program XPP show that a stable state in a chain of 
four oscillators coupled with d = 1.5 occurs close to x = (xi, £2, £3, £4) = (-6.4114408, 
7.3696656, 8.1984129, 7.3696656). Notice that the orbit of @ describing this state is not 
symmetric in the sense of proposition 3.4 although it is mapped into itself under R. This 



numerically obtained solution is used as an initial guess in Mathematica's implementation 
of Newton's method to determine the roots x^ of the set of equations 



a(x l 



(3rxi 



■) — 1.5(x( i _ 1 ) mo( j 4+1 — 2Xi + mod 4+l) 



1,2,3,4. (38) 



14 



Although Mathematica can be instructed to find roots to an arbitrary precision, so far only 
floating point arithmetic was used so that none of the results are rigorous yet. 

At this point all the quantities in the computations are redefined as intervals rather than 



floating point numbers and using interval arithmetic and theorem 6.2 we find an interval 
around each x l a in which the roots of equation ( |38| ) must lie. These bounds are now rigorous. 

In section [| it was shown that a fixed state and all of its X-multiples are stable if the 
polynomial px> (A, t) corresponding to that state has roots with negative real part for all 
t E [0, 7r/2]. Using interval arithmetic it is shown that 

pv(X,t) =3.61524 x 10 12 - 5.02753 x 10 7 cos At + (9.49583 x 10 11 - 1.33199 x 10 7 cos4f) A + 
(2.67802 x 10 11 - 4.94338 x 10 6 cos4i) A 2 + (4.90382 x 10 10 - 771725 cos At) A 3 + 



5 



(7.551 x 10 J - 144880cos4t) A + (9.81388 x 10 s - 13673. 3cos4t)A 
(1.03403 x 10 8 - 1560.66 cos 4i)A 6 + (9.39525 x 10 6 - 74~25cos4i)A 7 + 
(703005 - 5.0625 cos At)X 8 + 42026.7 A 9 + 2023.92 A 10 + 66.6667A 11 + A 12 

(39) 

for the fixed state under consideration. 

In the next step the roots of the polynomial px>(A, f) are found. To employ Mathemat- 
ical implementation of Newton's method we need a polynomial with coefficients that are 
floating point numbers rather than intervals. Floating point numbers inside the intervals 
which determine the coefficients of the polynomial pv(\ f ) can be chosen to define a poly- 
nomial pf>(A, |) which approximates Px>(A, §). The roots of {rf}}^ of Pp(A, are now 
found using Newton's method. 

The complex intervals^ containing the roots of pr>(A, t,) are found by using theorem 6.1 



We set zo equal to complex intervals around rf and use complex interval arithmetic to find 
the radius Kq and check the conditions of the theorem for each root r" . The 4 real roots 
and 4 complex pairs obtained by this procedure are given in the table below. 



-18.481 


-0.2463 ±i9.769 


-17.173 


-0.1961 ±i9. 317 


-15.437 


-0.1345 ±i9. 164 


-14.238 


-0.09054 ±i8. 622 



Since interval arithmetic was used in these calculations, these estimates are now rigorous. 
The remainder of the argument shows that these roots will not cross the imaginary axis as 
t varies. 

By lemma 5.1 the characteristic polynomial pv{\ t) takes the following form when eval- 
uated on the imaginary axis: 

12 11 

px>(ifJ.,t) = ^ (ay + /3j cos At) (i j (-l)z + i V] (atj + (3j cosAt)^ (-l) 1 ^ = 

j=0,j even j = l.j odd 



2 There are several ways in which complex intervals can be defined. In this case we use regions of the 
form z = x + iy where x and y are real intervals. These are simply rectangular regions in C. 
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Since the roots of px>(A, t) depend continuously on the parameter t, if px>(X,t) has roots 
with positive real part for some t\ then there must exist a to such that p-p(A, to) has a root 
on the imaginary axis. In other words there must exist a /iq such that 

p%(fx o ,t o )=p I v ( fM) ,to) = 0. (40) 

We will use interval arithmetic to show that this cannot happen. The polynomial px>(X,t) 
is split into a real and complex part to avoid using complex interval arithmetic in the 
numerical calculations since complex interval arithmetic leads to much rougher estimates 
than real interval arithmetic. 

By Gresgorin's theorem the eigenvalues of the matrix V(t) lie in the union of discs Cj = 
{z : \z — T>u\ < Ri} where T>ij are the entries in the matrix V(t) and Ri = J2j=i j^i l^ul- 
A direct computation shows that the intersection of Uj^C^ with the imaginary axis is 
contained in the interval [— 17i, 17i] and it is therefore sufficient to show that p%([i,to) and 
Px>(Mi *o) are not zero simultaneously for any value of t and fi G [— 17, 17]. 

Since the coefficients of Px>(A, t) are ? periodic in t, it is sufficient to show that (40) 
is not satisfied for any A G [—17,17] and t G [0, §]. This is shown by subdividing these 
intervals into a sufficient number of subintervals Jx n and i m so that [—17, 17] C uf =1 fi n and 
[0, ^] C U^ =1 i m with the property that the intervals p%(ft n ,tm) and p\j(p.n,t m ) do not 
contain zero simultaneously for any given pair of subintervals p, n and t m . Therefore the 
roots of pu(A, t) stay in the left half plane for all t G K and all i^-multiples of the fixed state 
under consideration are also stable. 

The paths of the roots of Pt>(\ t) as t is varied are shown in figure |ll|. 

The Mathematica code used in these calculations is available at 

http://math. bu. edu/ 'people /josic/ 'research/ code. 

Acknowledgment: The authors research was supported in part by NSF grant DMS- 
9803164. K.J. would like to thank the Department of Mathematics at Boston University for 
its hospitality while he conducted his research. The authors also thank J. P. Eckmann for 
useful discussions concerning this problem. 



A Proof of Lemma |b.l| 



Sublemma A.l Assume A is an (n + m) x (n + m) matrix of the form 

B D+ ... 01 

D- 

C 







(41) 



where m > n, B and C are arbitrary n x n, respectively m x m matrices and D are n x n 
matrices of the form 



D ± = 



d± ... 0" 
... 

... 



(42) 
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Writing the determinant of A as 



det A — sgn a J^J < 



then if a^ a {\) = oi, n +i = d + 



then 



M+n 



:,<7(i) 



(43) 



(44) 



onfo/ i/ a. 



n+l,<r(n+l) 



*n+l,l 



This lemma simply states that the only nonzero products in ( |43| ) containing d+ as a factor 
necessarily contain d_ as a factor. Notice that since we can interchange rows and columns 
when calculating determinants the opposite is also true: nonzero products containing d_ as 
a factor necessarily contain d + as a factor. 

Proof: Write the matrix A as 



A = 





i>l,2 


• • ■ &l,n d + 





&2,1 




• • • &2,„ 





d- 



























C 














(45) 



If 

a n+i,a(n+i) — a n+i,i = d+ and a(i) = 1 for some 2 < i < n another n— 1 nonzero factors 
in the product (Q) need to be chosen from columns 2 through n. However, barring rows 
1 and i which were already used, there are only n — 2 rows remaining which have nonzero 
entries in these columns. This shows that factors a,2,a(2)i a 3,o-(3)? ■ ■ ■ , &n,cr{ri) must come from 
columns 2 through n. Since the rows 1 through n are now exhausted, the only nonzero 
element remaining in column 1 is d_. This proves the lemma. 

Sublemma A. 2 The characteristic polynomial of the Km x Km matrix 



A = 



Mi E(t)+ 
E(t)~ M 2 E{t)+ 



E{t)~ M^-i E(t)+ 
E(t)- M K 



where Mj are arbitrary m x m matrices and E^t)^ are m X m matrices defined as 



de ±lt ... 0' 



does not depend on t. 
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Proof: Define the matrix C = A — XI. We need to show that the determinant of C is 
independent of t. The proof is by induction on K . When K = 1 the matrix C is itself f 
independent so the statement is trivial. The cofactors of the entries cij for 1 < j < m are 
of the form 



C2,l 


• ■ ■ C 2j 


■ ■ ■ Cl,m 




■ ■ ■ Cm,j 


■ ■ ■ C??i,m 


de~ u 


o ... 
































M 



(46) 



where the hat signifies the omission of that particular column. Since the columns and rows 



can be interchanged when the determinant is computed lemma A.l implies that any product 
in ( |43| ) containing de~ lt as a factor also contains ci jrl +i = as a factor and is thus zero. 
Lemma |A.l| also implies that any nonzero product containing ci in +i = de lt as a factor 
necessarily contains ci^+i = de~ %t as a factor and so does not depend on t either. Since 
this exhausts all possible cofactors of the nonzero elements of the first row the lemma is 
proved. <0> 

Lemma 5T can now be proved. We need to show that the determinant of C = B — XI 
is of the form (|35|). Definition ([||) will again be used for the determinant. We now have 
three cases: 

Case 1: Assume that for a in the product in ( [43|) we have 1 < <r(l) < m. Since the next 
m — 1 rows have zero entries in all but the first m columns this particular product will be 
zero unless it contains Cj t \ as a factor for some 1 < j < m. This excludes the possibility that 
either Civx-iw+i = der tl or ctK-i)m+i,i = de lt appear as factors. If we fix cr(l), . . . , c(m), 
the part of the sum ( [f3|) corresponding to all such permutations a is equal to 



£-1,(t(1) ■ • ■ Cm,£r(m) 



M 2 E{t)+ 
E{t)- M 3 E(t)+ 







15(f)- 




M K -i 
E{t)~ 






E{ty 

M K 



by cofactor expansion. By sublemma A. 2 these summands will not depend on t. 

Case 2: Assume that ct(1) = m + 1 so that the first factor in the product in 
ci.m+i = de li . Then the elements C2 i<T (2), ■ ■ • , Cm,a-(m) must come from the matrix 



(47) 



3b is 





C2,2 • 


C 2 , m 


c = 


C3,2 • 


C3,m 




Cm, 2 





by an argument identical to the one given in the second paragraph of sublemma (A.l) 
The only nonzero factors that can be chosen from column 1 are therefore c m+1 M (m+i) - 
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Cm+i : i = de lt and C(_ Kr _i) TO+ i )0 .(_ K: _i) ro+1 = C(x-i) m +i,i = de lt . We therefore have the 
following two subcases: 

Case 2a: If c m +\^i = de~ lt is chosen it cancels the term ci. m +i = de lt that was chosen 



as the first factor in the product. In this case we are again in the situation of sublemma A. 2 
using Laplace expansion as in case 1, and so this product will not depend on t. 

The contributions to the characteristic polynomial from Case 1 and 2a are time indepen- 
dent and therefore enter the term ^Sjf-iw+i o^A 1 in (|3^). The remainder of the lemma 
follows from induction. 

Case 2b: If the assumptions of Case 2 hold and C(K-i)m+i cr((K-i)m+i) = c (K-i)m+i,i = 
de lt then if a permutation a leads to a nonzero product in ( [Si) it must satisfy cr{jm + 1) = 
(j + l)m + 1 for all j < K — 1. This simply means that any nonzero product of this kind 
contains all the entries c jm+ i.( :)+ i) m+ i = de lt as factors. This assertion is proved as follows: 

The factor ci >m +i is contained in the product by assumption. Assume that a(lm + 1) = 
(I + l)m + 1 for all I < j. It will be shown that in this case a(jm + 1) = (jm + l)m + 1 
which will complete the induction argument. 

By assumption the cx((j — l)m + 1) = jm + 1. The only nonzero elements in columns 
jm + 2, . . . {j + 1)™ come from the submatrix 

Cjm+2,jm+2 ■ ■ ■ Cjm+2 : (j+l)m 
Cjm+3,jm+2 • ■ ■ Cjm+3,(j-t-l)m 



p(j+l)m,jm+2 • • • C(j + l)m,(j + l)m_ 



and since an entry from row jm + 1 has already been chosen, the coefficients Cj m+ 2,a(jm+2) > 
. . . , Cy + x)m,cr((i+i)m) must all come from this submatrix. The only nonzero entries re- 
maining on row jm + 1 are c jm+1Aj _ 1)m+1 = de~ lt and c jm+1 ^ j+1)m+1 = de lt . By the 
induction hypothesis we can exclude the first possibility and we are left to conclude that 
c jm+i.(j+i)m+i — de lt for any 1 < j < K is a coefficient in any nonzero product in the case 
2b. 

This shows that any permutation a of the type described in Case 2b contributes a factor 
r y (y e K ' lt r{\) to the characteristic polynomial (|35|). 

Case 3: The only remaining nonzero entry in row 1 of matrix C is c 1 (x-i)m+i — de~ lt . 
An argument identical to the one presented in case 2 shows that any permutation a leading to 
a nonzero summand in ( |43| ) must satisfy either C(K-i)m+i,o-((K-i)m+i) — C(K~i)in+i,i = de lt 
or C( i{ -_ 1 ) m+lj(T ((K-i)m+i) = C(K-i)m+i,(K-2)m+i = de~ lt which again leads to two subcases: 
Case 3a: If c^-i) m +i.cr{(K-i)m+ i) = de lt the situation is virtually identical to case 



2a. Laplace expansion and sublemma A. 2 show that the contributions to the characteristic 
polynomial are independent of t. 

Case 3b: If C(K~i)m+i.a((K-i)m+i) = de~ %t the situation is similar to case 2b. A parallel 
argument shows that all the terms de~ lt in the matrix must be coefficients in any nonzero 
product. Moreover due to the properties of the matrix, to each product ^ <y e Klt r{\) from 
case 2b corresponds a product r ) a e^ Klt r{\) from this case. Since these are the only time 
dependent contributions to the characteristic polynomial the lemma is proved. <C> 



References 



19 



[I] H.D.I. Abarbancl, et. al. "Synchronized Action of Synaptically Coupled Chaotic Neu- 
rons," Neural Comp. 8, 1567 - 1602 (1996). 

[2] V.S. Afraimovich, S-N Chow, Topological Spatial Chaos and Homoclinic Points of Z d - 
Actions in Lattice Dynamical Systems, Japan J. of Ind. and Appl. Math 12 367-383 
(1995). 

[3] V.S. Afraimovich, S-N Chow, Hyperbolic Homoclinic Point of Z d -Actions in Lattice 
Dynamical Systems, Int. J. of Bif. and Chaos 6 1059 - 1075 (1996). 

[4] Afraimovich, V. S., et. al. Stability, Structures and Chaos in Nonlinear Synchronization 
Networks, World Scientific (1994). 

[5] V.S. Afraimovich and W-W Lin, Synchronization in Lattices of Coupled Oscillators with 
Neumann/Periodic Boundary Conditions, preprint. 

[6] D.G. Aronson, G.B. Ermentrout and N. Kopell, Amplitude Response of Coupled Oscil- 
lators, Physica D 41 403-449 (1990). 

[7] PAshwin, J. Buescu, and I. Stewart, "From Attractor to Chaotic Saddle: a Tale of 
Transverse Instability," Nonlinearity 9, 703 - 737 (1996). 

[8] L. Bunimovich, "Localized Solutions in Lattice Systems and Their Bifucations Caused 
by Spatial Interactions," Nonlinearity 11 1539-1545 (1998). 

[9] L. Bunimovich, "Coupled map lattices: one step forward and two steps back," Physica 
D 86 248-255 (1995). 

[10] L. Bunimovich and Ya.G. Sinai, "Space-time Chaos in Coupled Map Lattices," Non- 
linearity 1 491-516 (1988). 

[II] Carvalho, A.N., Rodrigucs, H.M. and Dlotko, T., "Upper Semicontinuity of Attractors 
and Synchronization," to appear in J. Math. Anal, and App. 

[12] R.L. Devaney, "Homoclinic Bifurcations and the Area-Conserving Hcnon Mapping," J. 
Diff. Eq. 51 254-266 (1984). 

[13] G.B. Ermentrout, N. Kopell, "Oscillator Death in Systems of Coupled Neural Oscilla- 
tors," SIAM J. Appl. Math. 50 125-146 (1990). 

[14] E. Fermi, J. Pasta, and S. Ulam, "Studies of Nonlinear Problems," Los Alamos Report 
LA-1940. Reprinted in Lcc. App. Math 15, 143 - 156 (1955). 

[15] J.K. Hale, Diffusive Coupling, Dissipation and Synchronization, J. Dyn. and Diff. Eq. 
9 1 - 52 (1997). 

[16] J.F. Heagy, T.L. Carroll and L.M. Pecora, "Synchronous Chaos in Coupled Oscillator 
Systems," Phys. Rev. E 50, 1874 - 1885 (1994). 

[17] R.A. Horn and C.R. Johnson, Matrix Analysis, Cambridge University Press (1985). 



20 



[18] R. Huerta, M. Bazhenov, and M.I. Rabinovich, "Clusters of Synchronization and Insta- 
bility in Lattices of Chaotic Neurons," Europhy. Lett. 43, 719-724 (1998). 

[19] K. Josic Synchronization of Chaotic Systems, Ph.D. Thesis, Pennsylvania State Uni- 
versity (1999). 

[20] J. Moser, Proof of a Generalized Form of a Fixed Point Theorem Due to CD. Birkhoff 

[21] A.M. Ostrowski, Solutions of Equations and Systems of Equations, Academic Press 
(1960). 

[22] D.R. Orcndovici, Ya. B. Pesin, Chaos in Traveling Waves of Lattice Systems of Un- 
bounded Media, preprint. 

[23] M.I. Rabinovich, private communication (1999). 



21 



List of Figures 



1 A stable state of two breathers in a ring of 16 Lorenz systems coupled through 



the x variable at d x — 15. The first series of figures shows a snapshot of the 



variables of the different systems at different points of one oscillation. In 



these figures, the position of the oscillator in the ring is given on the .x-axis 



while Xj(t) is given on the y-axis. The timeseries of x\(t) and x^jt), which 



are nearly periodic, are given below 



23 



2 Typical timeseries of xi^jt) at coupling strengths d, x = 5, 15, 20, 25 from top 



to bottom. Many features of the system remain similar at different coupling 



values. The oscillations have a smooth envelope, a feature that is typical to 



solutions for coupling strengths at which the network behaves coherently, but 



is not synchronized exactly. Even the transients that precede a stable fixed 



state have similar shapes 24 



The system is very sensitive to initial conditions. The eventual stable state 



in a chain of 32 oscillators with d x = 40 is shown. Only one initial condition 



xq(0) is different for each figure - it was incremented by 0.01. Very different 



states arc reached. This suggests that the basins of attraction arc intertwined 



in a complicated way. The stable states in these figures are typical, although 



fixed states with one "hump" are more common than states with more "humps" 



25 



The wave when d x = dy_ = 6. This seems to be a globally stable solution ol 



the equations 26 





5 


The action of the two systems in equation (9) 






6 


The invariant circles ar 


Dund (0,0) give way to complicated behavio 


t close to 




(1/2, 1/2) due to the transversal intersection of W u (2/3, 2/3) and W (2/3, 2/3). 




7 


The region H is bounded by the two cubics F(D) and F(G) and D. The 






figure on the left W corresponds to the case d = 20. For d > 20 


f e these 






curves do not bound a compact region 




8 


If W u (p) coincides with the stable manifold of the origin then vi 


must be 
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Figure 1: A stable state of two breathers in a ring of 16 Lorenz systems coupled through 
the x variable at d x — 15. The first series of figures shows a snapshot of the x variables 
of the different systems at different points of one oscillation. In these figures, the position 
of the oscillator in the ring is given on the x-axis, while Xi(t) is given on the y-axis. The 
timeseries of x\(t) and x$(t), which are nearly periodic, are given below. 
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Figure 2: Typical timeseries of x 15(f) at coupling strengths d x = 5,15,20,25 from top to 
bottom. Many features of the system remain similar at different coupling values. The oscil- 
lations have a smooth envelope, a feature that is typical to solutions for coupling strengths at 
which the network behaves coherently, but is not synchronized exactly. Even the transients 
that precede a stable fixed state have similar shapes. 
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Figure 3: The system is very sensitive to initial conditions. The eventual stable state in a 
chain of 32 oscillators with d x = 40 is shown. Only one initial condition xo(0) is different 
for each figure - it was incremented by 0.01. Very different states are reached. This suggests 
that the basins of attraction are intertwined in a complicated way. The stable states in these 
figures are typical, although fixed states with one "hump" are more common than states 
with more "humps" . 
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Figure 4: The wave when d x = d y = 6. This seems to be a globally stable solution of the 
equations. 
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ure 5: The action of the two systems in equation 
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x= 1/2 



Figure 6: The invariant circles around (0,0) give way to complicated behavior close to (1/2, 
1/2) due to the transversal intersection of W"(2/3, 2/3) and VF s (2/3, 2/3). 
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Figure 7: The region if is bounded by the two cubics F(D) and F{G) and D. The figure 
on the left W corresponds to the case d — 20. For d > 20 + e these curves do not bound a 
compact region. 
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Figure 8: If W u (p) coincides with the stable manifold of the origin then vi must be tangent 
to W u (p) at the origin. Since v\ points outside of H this is impossible. 




Figure 9: A numerical approximation of the manifold W u {p) is shown on the right. The 
complex structure of the manifold suggests complicated behavior on some invariant set of 
points. 
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Figure 10: The fixed state that will be used as an example, and its 4-multiple. 
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Figure 11: A numerical computation of the paths that the 8 complex eigenvalues with 
smallest real part trace out as t varies. Only one eigenvalue in the complex conjugate pair 
is shown. Notice the small trace of the eigenvalue with the least negative real part. 
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